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mance of a quantum computer, it is necessary to be able to perform operations 
fault-tolerantly on encoded states. I present a general theory of fault-tolerant 
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operations based on symmetries of the code stabilizer. This allows a straight- 
forward determination of which operations can be performed fault-tolerantly 
on a given code. I demonstrate that fault-tolerant universal computation is 
possible for any stabilizer code. I discuss a number of examples in more detail, 



including the five-qubit code. 
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The development of quantum error-correcting codes [T|-§[ has stirred great hopes for 
conquering errors and decoherence in quantum computers. However, just the existence of 
codes, even very good codes, is not sufficient. It is also necessary to be able to perform 
operations on encoded states without a catastrophic spread of existing errors. However, 
until now, it was only known how to implement a universal set of gates fault-tolerantly on 
a few codes 

While a general formalism for quantum error-correcting codes is known, the code stabi- 
lizer [£§-|6[], the theory of fault-tolerant operations is in a great deal of disarray. A quantum 
gate, unlike a classical gate, can cause errors to spread both forwards and backwards through 
the gate. The goal of fault-tolerant operations is to prevent the spread of errors within a 
block, which could change a single correctable error into two errors, which is perhaps more 
than the code could handle. Even if we use codes that correct more than one error, the 
spread of errors within a block rapidly reduces the code's tolerance for errors. Therefore, 
I define a fault-tolerant operation to be one for which a single operational error can only 
produce one error within a single encoded block. The assumption is that storage errors on 
different qubits are independent and that gate errors can only affect qubits which interact 
via that gate. 

A transversal operation, in which the operation acts independently on each qubit in the 
block, is a prototypical fault-tolerant operation. For instance, a bitwise controlled NOT 
operation from one block to another is fault-tolerant, since errors can spread only between 
corresponding qubits in the two blocks. 

Unfortunately, most bitwise operations applied to most codes will not map one valid 
codeword to another. Until now, there has been no general theory of what a bitwise operation 
will do to the codewords of a given code. I show below that a bitwise operation will transform 
the stabilizer of a code. If the stabilizer is rearranged, but otherwise left unchanged, the 
operation will take codewords to codewords. This will give us a few basic operations on 
various codes with which to start our analysis. 

In the quest to perform universal quantum computation, we are not limited to unitary 



operations. We can also perform measurements. In section |HT[ I analyze the behavior of 
certain states when a measurement is made. This allows us to see what operations we can 
derive from the basic operations by using ancillas and making partial measurements of the 
state. Ultimately, this will allow us to perform universal computation on a large set of codes, 
among them the five-qubit code, a class of distance 2 codes, and the code encoding 3 qubits 
in 8 qubits. 

I. ENCODED NOT AND PHASE 

Before I advance into the full theory of fault-tolerant operations, I will discuss how to 
perform encoded NOT and Phase gates on any stabilizer code. The behavior of these gates 
under more general transformations will tell us what those transformations actually do to 
the encoded states. 

The stabilizer S of a code is an abelian subgroup of the group Q generated by the 

operations 
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acting on each of the n qubits of the code. I will sometimes write Q n to explicitly indicate 
the number of qubits acted on. The codewords of the code are the states \ip) which are 
left fixed by every operator in S. Operators in Q which anticommute with some operator 
in S will take codewords from the coding space into some orthogonal space. By making 
a measurement to distinguish the various orthogonal spaces, we can then determine what 
error has occured and correct it. A quantum code encoding k qubits in n qubits will have a 
stabilizer with n — k generators. 

However, there are, in general, a number of operators in Q that commute with all of the 
operators in S. The set of such operators is the normalizer N(S) of S in (/.[] S is itself 



1 Strictly speaking, this is the centralizer of S, but in this case it is equal to the normalizer, since 



contained in the normalizer, but in general the normalizer is larger than just S. If S contains 
2 a operators (so it has a generators), the normalizer will be generated by 2n — a = n + k 
operators. (Q has a total of 2 2n elements, half of which will commute with any other fixed 
element of Q.) 

The elements of the normalizer will change one codeword to another, and therefore have a 
natural interpretation in terms of encoded operations on the code words. Suppose we extend 
the stabilizer into a maximal set of n commuting operators. Then consider those codewords 
which, besides being +l-eigenvectors of the stabilizer generators, are also eigenvectors of 
the additional k operators from N(S). Let these codewords be the basis codewords for our 
code, giving the encoded |0 . . . 00), |0 . . . 01), ... |1 ... 11). The state which has eigenvalue 
+ 1 for the k new operators will be the encoded |0 . . . 0), the state which has eigenvalue —1 
for everything will be the encoded |1 . . . 1), and so on. Then the k new operators have the 
interpretation of being the encoded Z operators on the k encoded qubits. We will write 
these encoded Z operators as Z t for the ith encoded qubit, or Z when there is just one 
encoded qubit. 

Now, the remaining elements of N(S) will not commute with all of the encoded Z 
operators. We can choose a set of generators for N(S) such that each of the last k operators 
commutes with everything except a single one of the Z operators. These generators are 
then the encoded bit flip operators X; t (or X when there is just one). An arbitrary element 
of N(S) is some other encoded operation. If two elements of N(S) differ by an element of 
the stabilizer, they act the same way on any code word (since the stabilizer element just 
fixes the codeword). Therefore the actual set of encoded operations represented in N(S) is 
N{S)/S. 

DiVicenzo and Shor showed how to perform syndrome measurement and error correction 
fault-tolerantly on any stabilizer code pi. Using the same methods, we can measure the 
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eigenvalue of any operator in Q, even if it is not in S. This also enables us to prepare the 
encoded zero state of any stabilizer code by performing error correction and measuring the 
eigenvalue of the Z operators. 

II. MORE GENERAL OPERATIONS 

So far, we have only considered applying products of X, Y, and Z to the codewords. 
However, this is not the most general thing we could do. Suppose we have some totally 
arbitrary unitary transformation U we wish to apply to our codewords. How does this affect 
other operators, such as the elements of S and N(S)1 

UM\ip) = UMU^U\ip} 7 (2) 

so \ip) is an eigenvector of M if and only if U\ip) is an eigenvector of UMU^ . Furthermore, 
they have the same eigenvalue. Thus, by applying U to \ip), we effectively transform any 
operator M of interest into UMU^ . In order for the state \ip) to remain a codeword, the state 
U\ip) must still be in the coding space, so UMW must also fix all the codewords \ip) for every 
M e S. Let us consider a restricted set of possible Us, those for which UMW is actually 
in Q (so U is in the normalizer N(Q) of Q in U(n)). We will see that N(Q) is generated by 
Hadamard rotations, tt/2 phase rotations, and controlled NOT operations ||5|JTT|. Then, by 
the definition of the stabilizer and the coding space, we need UMU^ to actually be in S for 
all M G S. Therefore, U is actually in the normalizer of S in U(n). The same criterion was 
found previously by Knill ]ll|. Note that the normalizer of S in U(ri) is not necessarily a 
subset of N(G). 

When we restrict our attention to operations that are in both the normalizer of Q in 
U(n) and the normalizer of S in U(ri), it becomes straightforward to determine the oper- 
ation actually performed on the encoded states. First, note that the X and Z operators 
transform into operators that also commute with everything in S. Thus, we can rewrite 
them as products of the original Xs, Zs, and elements of S. The elements of S just give 



us the equivalence between elements of N(S) discussed in section |, so we have deduced 
a transformation of the encoded X and Z operators. Furthermore, we know this encoded 
transformation also lies in the normalizer of Qk- 

Typically, we want to consider transversal operations U, which are equal to the tensor 
product of single-qubit operations (or operations that only affect one qubit per block). For 
the moment, we will only consider operations of this form and see what collections of them 
will do to the stabilizer. Before launching into an analysis of which gates can be used on 
which codes, I will present an overview of the gates that are amenable to this sort of analysis. 

For instance, one of the simplest and most common fault-tolerant operations is the 
Hadamard rotation 

fl= ^C -J- (3) 

Let us see what this does to X, Y, and Z. 

^C -i)(-i !) = Q = * (5) 

«™h(; _;)(■; ;)-(_■ ;)--* <«> 

Therefore, applying R bitwise will switch all the Xs and all the Zs, and give a factor of 
—1 for each Y . If we do this to the elements of the stabilizer and get other elements of the 
stabilizer, this is a valid fault-tolerant operation. The seven-qubit code is an example of a 
code for which this is true. 

Another common bitwise operation is the % phase 

1 0\ 

(7) 
if 

On the basic operations X, Y, and Z it acts as follows: 

/l 0\ /0 -i\ /0 -f 

PXF, io J(, o)i, o 1 ^ (8) 
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This switches X and Y", but with extra factors of i, so there must be a multiple of 4 As 
and Ys for this to be a valid operation. Again, the seven-qubit code is an example of one 
where it is. Note that a factor of i appears generically in any operation that switches Y 
with X or Z, because Y 2 = — 1, while X 2 = Z 2 = +1. The operations in N(Q) actually 
permute ax = X, az = Z, and ay = iY, but for consistency with earlier publications I have 
retained the notation of X, Y, and Z. The most general single qubit operation in N{Q) can 
be viewed as a rotation of the Bloch sphere permuting the three coordinate axes. 

We can also consider two-qubit operations, such as the controlled NOT. Now we must 
consider transformations of the two involved blocks combined. The stabilizer group of the 
two blocks is S x S, and we must see how the basic operations X <g) I, Z ® I, I ® X, 
and I <S> Z transform under the proposed operation. In fact, we will also need to know the 
transformation of X <g> Y and other such operators, but the transformation induced on Q x Q 
is a group homomorphism, so we can determine the images of everything from the images 
of the four elements listed above. 

It is straightforward to show that the controlled NOT induces the following transforma- 
tion: 

X <g> i" -► X ® X 

Z®I^Z®I (11) 

I®X -f I®X 

I®Z^Z®Z. 

It is easy to see here how amplitudes are copied forwards and phases are copied backwards. 
The transformation laws for R, P, and CNOT are also given in ||. 

There are a number of basic gates in N(Q) beyond the ones given above. As with the 



examples above, any gate can be characterized by its transformation of the generators of Q 
(or Q x Q for two-qubit operations, and so on). The primary constraint that must be met is 
to preserve the algebraic properties of the operators. In fact, there is a complete equivalence 
between the possible gates and the automorphisms of Da, (the group of products of /, X, Y, 
and Z) or direct products of copies of D 4 (for multiple-qubit gates) [|TB|]. 

Given any such automorphism, we first substitute iY for Y to get the actual transforma- 
tion. Then we note that |0) is the "encoded zero" for the "code" with stabilizer {/, Z}. We 
know how Z transforms under U, so |0) transforms to the state fixed by UZU^ . In addition, 
|1) = -A|0), so U\l) = UXWU\0). For instance, consider the cyclic transformation 

T = X^iY^Z^X. (12) 

Since Z — > X, 

|0)-l/v / 2(|0) + |l)). (13) 

Also, X —* iY, so 

|1> -, i/V2 Y(\0) + |1» = -i/y/2 (|0) - |1». (14) 

Thus, the matrix for T is 

T ^(! ")■ (15) 

We can perform a similar procedure to determine the matrix corresponding to a multiple- 
qubit transformation. 

The next question of interest is how much have we restricted our computational power 
by restricting our attention to the normalizer of Q? Again, the normalizer of Q is exactly 
the group generated by the Hadamard transform R, the phase P, and the controlled-NOT. 
I will prove this in section ^TJ. Unfortunately, this group alone is of only limited interest. 
Knill |I3| has shown that a quantum computer using only operations from this group can 



be simulated efficiently on a classical computer .0 However, the addition of just the Toffoli 
gate to this group is sufficient to make the group universal [0]. 

III. MEASUREMENTS 

Now I will discuss what happens if we perform a measurement on a stabilizer code. 
Measuring individual qubits of an actual code is not of great interest, but the results of this 
section will be quite helpful in determining what can be done by combining measurements 
and specific fault-tolerant operations. 



Now, using the method of DiVincenzo and Shor [ 10| , we can measure any operator A in 
Q. There are three possible relationships between A and 5*. First of all, A could actually be 
in S. Then measuring A tells us nothing about the state of the system and does not change 
it at all. The result of this measurement will always be +1 for a valid codeword. The second 
possibility is for A to commute with everything in S but not to actually be in S. Then A 
is equivalent to a nontrivial element of N(S)/S and measuring it will give us information 
about the state of the system. This is usually inadvisable. 

The third possibility, that A anticommutes with something in S, is the most interesting. 
In this case, we can choose the generators of S so that A anticommutes with the first 
generator Mi and commutes with the remaining generators M 2 , . . . , M n _k (we can do this 
since if generator Mj anticommutes with A, we can replace it with MiMj, which commutes). 
Then measuring A does not disturb the eigenvectors of M 2 through M„_fe, so those still fix 



2 The argument goes as follows: we start with an n-qubit state |0) which is the single state for 
the stabilizer code (Z\, . . . , Z n ). Each operation transforms the state and the stabilizer as above. 
We can follow each transformation on a classical computer in 0(n 2 ) steps. A measurement picks 
at random one of the basis kets in the codeword, which can also be chosen classically ||,14|. This 



still leaves the question of partial measurement of the full state, but the results of section III show 
that this can also be classically simulated. 



the new state, and are in the new stabilizer. The eigenvectors of Mi are disturbed, however, 
and Mi no longer fixes the states. Measuring A applies one of the projection operators P + 
or P_, where 

P± = \(I±A). (16) 

Then M\P_M X = M\MiP + = P + , so if \if>) is some codeword, 

M{PJ$) = AfJP_Mi|^) = P+\i/j). (17) 

If the measurement result is +1, we do nothing else, and have thus applied P + . If the 
measurement result is — 1, apply M\ = Mi, resulting in the overall application of P + . Either 
way, A fixes the new state. This puts the system into the space with stabilizer generated 
by A, M 2 , . . . , M„_fc. From now on, I will often say "measure" when I mean "measure and 
correct for a result of — 1." 

Note that this construction works outside the framework of stabilizer codes. All we really 
need is a state \ip), with M\ip) = \ip) for some unitary M. Then, as above, we can perform 
the projection P + for any operator A satisfying A 2 = 1 and {M, A} = 0. Note that if A is 
some local measurement, either \ip) is not entangled between the region affected by A and 
the region unaffected by it, or M is a nonlocal operator. 

We will want to know just where in the space a given state goes. To do this, look at the 
elements of N(S)/S. If the state starts in an eigenvector of N, it will also be an eigenvector 
of N' = MN for all M G S. After measuring A, the state will no longer be an eigenvector 
of N if N anticommutes with A, but it will still be an eigenvector of MiN, which commutes 
with A. Furthermore, the eigenvalue of M\N stays the same. Therefore, by measuring 
A (and correcting the state if the result is —1), we effectively transform the operator N 
into MiN . We could equally well say it is transformed to MMiN instead, where M G S 
commutes with A, but this will produce the same transformation of the cosets of N(S)/S 
to N(S ; )/S' (where S' is the stabilizer after the measurement). Of course, if iV commutes 
with A, measuring A leaves iV unchanged. 

10 



Let us see how all this works with a simple, but very useful, example. Suppose we have 
two qubits, one in an arbitrary state \i/j), the other initialized to |0). The space of possible 
states then has stabilizer I (g) Z. Suppose we perform a controlled-NOT from the first qubit 
to the second. This transforms the stabilizer to Z ® Z. Now let us measure the operator 
I ® iY (we use the factor of i to ensure that the result is ±1). This anticommutes with 
Z <g) Z, so if we get +1, we leave the result alone, and if we get —1, we apply Z <g> Z to the 
state. The new state is in a +l-eigenstate of I <8> iY, that is, |0)(|O) + i\l)). 

How is \ip) related to \<f>)l For the original "code," X = X ® I and Z = Z ® I. After 
the CNOT, X = X <g> X and Z = Z <g) I. X <g> X does not commute with J <8> zF, but the 
equivalent operator Y ®Y = {X <g> X)(Z ® Z) does. Z ® I does commute with J ® iY, so 
that stays the same. Since the second qubit is guaranteed to be in the +1 eigenstate of iY, 
we might as well ignore it. The effective X and Z operators for the first qubit are thus — iY 



and Z respectively. This means we have transformed X — > —iXZ and Z — > Z. This is the 
operation PL 

This example is simple enough that it is easy to check: 

|00>-|00> = |0>i[(|0>+<|l» + (|0>-i|l»] (18) 

-|0>(|0)±i|l» (19) 

-|0)(|0)+i|l» (20) 

|10>-|ll> = |l)i[-(|0>+i|l)) + (|0>-»|l))] (21) 

-i|l>(=F|0>-i|l» (22) 

-±i|l>(=F|0>=F*|l» = -»11>(|0)+i|l)). (23) 

Thus, ignoring the second qubit gives |0) — > |0) and |1) — > — i\l), which is PL 

This result is already quite interesting when coupled with the observation that P and 
CNOT suffice to produce R as long as we can prepare and measure states in the basis |0) ± 1 1) 
||. To do this we start out with the state \ip) plus an ancilla |0) + |1). Thus, the initial 
stabilizer is I ® X, X = X <g> /, and Z = Z <g> I. Apply a CNOT from the second qubit to 
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the first. Now the stabilizer is X <g> X, X = X <g> /, and Z = Z ® Z . Apply P to the second 
qubit, so the stabilizer is X <g> iY, X — X <S> I, and Z = Z ® Z . Measure I ® X, performing 
X ® iY if the result is —1. This produces X = X ® I and Z = iY <g> X, so dropping the 
second qubit results in the transformation Q: X — *■ X, Z — *■ iY. But i? = PQ^P: 

X -»• «y -»• Z -»■ Z (24) 

Z ^ Z ^ -iY -► X (25) 

Coupled with the previous result, which derives P from CNOT, this allows us to get any 
single qubit transformation in the normalizer of Q provided we can perform a CNOT oper- 
ation. 

Another interesting application is to gain a new viewpoint on quantum teleportation. 
Suppose we have three qubits which start in the state |?/>)(|00) + |11)). The initial stabilizer 
is I <S> X <g) X and / <g> Z <g> Z, X = X <g> / <g> /, and ~Z = Z®I®I. We assume the third qubit 
is far away, so we can do no operations interacting it directly with the other two qubits. 
We can, however, perform operations on it conditioned on the result of measuring the other 
qubits. We begin by performing a CNOT from qubit one to two. The stabilizer is now 
I <g> X <g> X and Z<g)Z<g)Z, X = X<g)X<g)/, and Z = Z <g> I ® I. Measure X for qubit 
one and discard qubit one. If the measurement result was +1, we leave the state alone; if it 
was —1, we perform Z on qubits two and three. The stabilizer is now X®X, X = X£g>/ 
and Z = Z <g> Z . Now measure Z for the new first qubit. If the result is +1, we leave the 
final qubit alone; if it is — 1, we apply X to the last qubit. This results in X = X and 
Z = Z, both acting on the last qubit. We have succesfully teleported the state \ip). The 
operations conditioned on measurement results needed for teleportation arise here naturally 
as the corrections to the stabilizer for alternate measurement results. The formalism would 
have told us just as easily what operations were necessary if we had begun with a different 
Bell state or a more complicated entangled state (as long as it can still be described by a 
stabilizer). 

I claimed before that products of R, P, and CNOT actually gave us all of the elements 
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of N(Q), and I am now ready to prove that. The one-qubit operations correspond to the six 
automorphisms of -D4 given by R, P, Q, T, T 2 , and of course the identity. We have already 
seen that Q = P^RP^ . Also, T = PQ\ so all one-qubit operations are covered. 

We can also perform all two-qubit operations. Every automorphism of -D4 x D4 can be 
produced by a composition of controlled NOT and single-qubit operations. For instance, 
take 

Z®I -»• x ®x 

I®Z^Z®Z (26) 

X ® I -»• %Y ® X 

I®X ^iZ®Y. 

This permutation can be produced by performing the cyclic permutation X — ► iY — *■ Z — ► X 
on the first qubit and a phase rotation X — > iY on the second qubit, and then performing 
a standard controlled NOT from the first qubit to the second qubit. It is straightforward to 
consider the other possibilities and show that they too can be written using a CNOT and 
one-qubit gates. 

I will show that the larger gates can be made this way by induction on the number of 
qubits. Suppose we know this to be true for all n-qubit gates, and we have an (n + l)-qubit 
gate U. On an arbitrary input state |0)|^>) + |1)|0) (where \ip) and \<f>) are n-qubit states), 
the output state will be 

(|0>|^> + |1>|^» + (10)10!) + |1)|0 2 ». (27) 

Suppose that under the applied transformation, M = U(Z ® I ® ■ ■ ■ ® I)W anticommutes 
with Z®I®- ■ -®I. If it does not, we can apply a one-qubit transformation and/or rearrange 
qubits so that M — X ® M', where M' is an n-qubit operation. Suppose we apply U to 
|0) \ip). If we were then to measure Z for the first qubit, we would get either 0, in which case 
the other qubits are in state \ipi), or 1, in which case the remaining qubits are in state IV^)- 
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The above analysis of measurements shows that \ipi) and \ip2) are therefore related by the 
application of M' . Define U' by U'\ip) = \ipi). Then 

U(\0)\4>)) = (I + M)(\0)®U>\4>)). (28) 

Let A" = U(X ®I <g> ■ ■ -®I)U' . Again, we can apply a one-qubit operation so that either 
N = Z®N' or N = I <g)N'. We can always put M and A" in this form simultaneously. Then 

[7(|l)|0» = iV[/(|O>|0» (29) 

= AT(I + M)(|O)(g)[/ / |0)) (30) 

= (I-M)N(\0)®U'\(P)) (31) 

= (I-M)(\0)®N'U'\<f>)), (32) 

using the above form of A" and the fact that {M, N} = 0. 

Now, U' is an n-qubit operation, so we can build it out of R, P, and CNOT. To apply 
U, first apply U' to the last n qubits. Now apply N' to the last n qubits conditioned on the 
first qubit being 1. We can do this with just a series of CNOTs and one-qubit operations. 
Now apply a Hadamard transform to the first qubit. This puts the system in the state 

(|0> + |1» <g) U'\tl)) + (|0> - |1» ® N'U'\<j>). (33) 

Now, apply M' to the last n qubits conditioned on the first qubit. Again, we can do this 
with just CNOTs and one-qubit operations. This leaves the system in the state 

|0) ® U'\i)) + |1) ® M'U'\ip) + |0) ® N'U'\(j>) - |1) ® M'N'U'\(j)} (34) 

= |0) ® U'\ij)) + M(|0) ® U'\vl))) + |0) ® N'U'\(f)) - M(|0) ® iV'C/'|0)) ( 35 ) 

= (/ + M)(|0) ® C/'IV')) + (/ - M)(|0) ® iV'C/'|0)), ( 36 ) 
which we can recognize as the desired end state after applying U. 

IV. OPERATIONS ON CSS CODES 

In this section, I will finally begin to look at the problem of which gates can be applied 
to specific codes. One of the best classes of codes for fault-tolerant computation are the 
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TABLE I. The stabilizer and encoded X and Z for the seven-qubit code. 

Calderbank-Shor-Steane (CSS) codes [@|§, which are converted from certain classical codes. 
These codes have a stabilizer which can be written as the direct product of two sectors, 
one of which is formed purely from Xs and one formed just from Zs. These two sectors 
correspond to the two dual classical codes that go into the construction of the code. 

Shor showed that a punctured doubly-even self-dual CSS code could be used for 
universal computation. An example of such a code is the seven-qubit code, whose stabilizer 
is given in Table |. From the stabilizer, we can now understand why such codes allow the 
fault-tolerant implemention of the Hadamard rotation, the it/2 rotation, and the controlled 
NOT. 

The Hadamard rotation switches X and Z. For a CSS code, this is a symmetry of the 
stabilizer if and only if the X sector of the stabilizer is the same as the Z sector. Therefore 
the two classical codes had to be identical, and the quantum code must be derived from a 
classical code that contains its own dual. As we can see, this works for the seven-qubit code. 
In order to understand what the Hadamard rotation does to the encoded states, we must 
look at what it does to the encoded X and Z operations. For a punctured self-dual CSS 
code, the X and Z operations can again be taken to be the same, so the Hadamard rotation 
will just switch them. It is therefore an operation which switches encoded X with encoded 
Z, and is thus an encoded Hadamard rotation. 
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Similarly, for a self-dual code, the n/2 rotation will convert the X generators into the 
product of all Ys. This just converts the X generators into their product with the corre- 
sponding Z generator, so this is a valid fault-tolerant operation, provided the overall phase 
is correctly taken care of. There is a factor of % for each X, so there must be a multiple 
of 4 Xs in each element of the stabilizer for that to work out in general. This will only be 
true of a doubly-even CSS code, which gives us the other requirement for Shor's methods. 
Again, we can see that the seven-qubit code meets this requirement. Such a code will have 
3 mod 4 Xs in the X operation, so the bitwise 7r/4 converts X to —iY. This is thus an 
encoded — 7r/2 rotation. 

Finally, we get to the controlled NOT. This can be performed bitwise on any CSS code. 
We must look at its operation onM®/ and I®M. In the first case, if M is an X generator, 
it becomes M ® M. Since both the first and second blocks have the same stabilizer, this is 
an element of S x S. If M is a Z generator, M <S> I becomes M <S> I again. Similarly, if M is 
an X generator, I <g> M becomes I ® M, and if M is a Z generator, I <g> M becomes M <g> M, 
which is again in S x S. For an arbitrary CSS code, the Xi operators are formed from the 
product of all Xs and the Z$ operators are formed from the product of all Zs. Therefore, 

X~(g) I -h> X~<g>X~ 

~Zi®I ^%®I (37) 

I®~X~i—> I (g)X~ 

I ®~Z~^~Z~®~Z~. 

Thus, the bitwise CNOT produces an encoded CNOT for every encoded qubit in the block. 
In fact, we can now easily prove that codes of the general CSS form are the only codes 
for which bitwise CNOT is a valid fault-tolerant operation. Let us take a generic element 
of the stabilizer and write it as MN, where M is the product of Xs and N is the product 
of Zs. Then under bitwise CNOT, MN <g> I -»■ MN <g> M, which implies M itself is an 
element of the stabilizer. The stabilizer is a group, so N is also an element of the stabilizer. 
Therefore, the stabilizer breaks up into a sector made solely from Xs and one made solely 
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TABLE II. The stabilizer and encoded X and Z for the five-qubit code, 
from Zs, which means the code is of the CSS type. 



V. THE FIVE-QUBIT CODE 



One code of particular interest is the five-qubit code ||lT| , p~5f j which is the smallest possi- 
ble code to correct a single error. Until now, there were no known fault-tolerant operations 
that could be performed on this code except the simple encoded X and encoded Z. One pre- 
sentation H of the five-qubit code is given in Table 0. This presentation has the advantage 
of being cyclic, which simplifies somewhat the analysis below. 

This stabilizer is invariant under the transformation T : X —* iY — ► Z —>■ X bitwise. 
For instance, 

Mi = X ® Z ® Z ® X ® I ^ -Y ® X ® X ®Y ® I = M 3 Mi. (38) 

By the cyclic property of the code, M 2 through M 4 also get transformed into elements of 
the stabilizer, and this is a valid fault-tolerant operation. It transforms 

X -> iY ->• Z. (39) 

Therefore, this operation performed bitwise performs an encoded version of itself. Oper- 
ations which have this property are particularly useful because they are easy to apply to 
concatenated codes |§|],[16 . 
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There is no nontrivial two-qubit operation in the normalizer of Q that can be performed 
transversally on this code. However, there is a three-qubit transformation T 3 that leaves 
S x S x S invariant: 



X®I®I 

Z®I®I 
I®X ®I 
I®Z®I 



iX ®Y ® Z 
iZ ®X ®Y 
iY ®X ® Z 
iX® Z ®Y 



(40) 



I ® I ®X ^ X ®x ®x 
I®I®Z->Z®Z®Z. 

On operators of the form M ® I ® I or I ® M ® I, this transformation applies cyclic 
transformations as above to the other two slots. Operators I®I®M just become M®M®M, 
which is clearly in S x S x S. The matrix of T3 is (up to normalization) 
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(41) 



As with T, this operation performs itself on the encoded states. A possible network to 
produce this operation (based on the construction in section p| ) is given in figure [l|. 

If we add in the possibility of measurements, this three-qubit operation along with T 
will allow us to perform any operation in the normalizer of Q. I will describe how to do this 
on unencoded qubits, and since T and T3 bitwise just perform themselves, this will tell us 
how to do the same operations on the encoded qubits. 



-© — ? — © — j — r-^h-^- 
o CD — i — iOTH> 



■A- 



i?— 



FIG. 1. Network to perform the T3 gate. 

To perform P, first prepare two ancilla qubits in the state |00) and use the data qubit 
as the third qubit. The original stabilizer is Z <8> I <8> I and I <g> Z <S> / , X = / ® / <g> X, 
and Z = I ® I ® Z. Now apply T3, so that the stabilizer is zZ <g> X <g> K and zX ® Z (g) K, 
X = X £g> X ® X, and Z = Z ® Z ® Z. Measure Z for the second and third qubits. The 
resulting X = iY ® I ® Z and Z = Z <g> Z <g> Z. Dropping the last two qubits, we have 
X — > iY and Z ^ Z, which is P. Again, Q = T^P and R = PQ^P, so we can perform any 
single qubit operation. 

To get a two-qubit operation, prepare a third qubit in the state |0) and apply T 3 . This 
results in the stabilizer Z<g>Z<g>Z, ~X[ = iX ®Y ® Z , X^ = iY (g) X (g) Z, ~Z[ = iZ ® X ®Y , 
and Z 2 = iX ® Z ($>Y . Measure X for the second qubit and throw it out. This leaves the 
transformation 

X <g) I — ► iY ® I 

I®X ^iY ®Z (42) 

Z ® I ^ iZ ®Y 
I®Z^iY®X. 

This operation can be produced by applying Q to the second qubit (switching Z and iY), 
then a CNOT from the second qubit to the first one, then P to the first qubit and T 2 
to the second qubit. Therefore, we can also get a CNOT by performing this operation 
with the appropriate one-qubit operations. This allows us to perform any operation we 
desire in the normalizer of Q. Note that section provides us with another way to get 
these operations. Having two methods available broadens the choices for picking the most 
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efficient implementations. 

In order to perform universal computation on the five-qubit code, we must know how to 
perform a Toffoli gate. Shor ]7| gave a method for producing a Toffoli gate that relied on 
the ability to perform the gate 

\a)\b}\c}^(-ir^\a}\b}\c), (43) 

where \a) is either |0 . . . 0) or |1 . . . 1) and |6) and \c) are encoded 0s or Is. For the codes Shor 
considered, this gate could be performed by applying it bitwise, because the conditional sign 
could be applied bitwise. All of the qubits in the first block are either or 1, so a controlled 
conditional sign from the first block will produce a conditional sign on the second two blocks 
whenever the first block is 1. 

For the five-qubit code, this gate is not quite as straightforward, but is still not difficult. 
To perform the two-qubit conditional sign gate on the five-qubit code, we need to perform a 
series of one- and three-qubit gates and measurements. However, if we perform each of these 
gates and measurements conditional on the value of a, we have performed the conditional 
sign gate on \b)\c) if and only if the first block is 1. 

From this, the rest of Shor's construction of the Toffoli gate carries over straightforwardly. 
It involves a number of measurements and operations from the normalizer of Q. We have 
already discussed how to do all of those. The one remaining operation that is necessary is 

\a)\d) -► (-l) ad \a)\d), (44) 

where \d) is an encoded state and \a) is again all 0s or all Is. However, this is just Z applied 
to \d) conditioned on the value of a, which can be implemented with a single two-qubit gate 
on each qubit in the block. Therefore, we can perform universal fault-tolerant computation 
on the five-qubit code. 

Note that there was nothing particularly unique about the five-qubit code that made the 
construction of the Toffoli gate possible. The only property we needed was the ability to 
perform a conditional sign gate. 
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FIG. 2. Network to perform the four-qubit gate. 



VI. GATES FOR ANY STABILIZER CODE 



Consider the following transformation: 

X <g> / ® I ® I - 

I®X®/®I- 

J<g/<gX®/- 

J®/®/®x- 



X (g) X <g> X (g) / 

I (g) X (g) X (g) X 
X (g I (gX (gX 
X (gX (g I (gX 



(45) 



Z®I®I®I->Z<®Z®Z®I 

I(gZ(gI(gI^I(gZ(gZ(gZ 
I®I(gZ®I^Z(gI(gZ(gZ 
I(gI(gI®Z^Z(gZ(gI(gZ. 

A possible gate array to perform this operation is given in figure |2|. This operation takes 
M(gJ(g/(gJ to M(gM(gM(g J, and cyclic permutations of this, so if M G S, the image of 
these operations is certainly in S x S x S x S. This therefore is a valid transversal operation 
on any stabilizer code. The encoded operation it performs is just itself. There is a family 
of related operations for any even number of qubits (the two-qubit case is trivial), but we 
only need to concern ourselves with the four-qubit operation. 

Suppose we have two data qubits. Prepare the third and fourth qubits in the state |00), 
apply the above transformation, and then measure X for the third and fourth qubits. The 
resulting transformation on the first two qubits is then: 
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I®X ^ / ® A (46) 

Z®I -»• Z®I 

i®z ^ z®z. 



This is precisely the controlled NOT. Since I showed in section [I I J that the CNOT was 
sufficient to get any operation in N(G), we can get any such operation for any stabilizer 
code! In fact, using the Toffoli gate construction from section [V], we can perform universal 
computation. 

Actually, this only gives universal computation for codes encoding a single qubit in a 
block, since if a block encodes multiple qubits, this operation performs the CNOT between 
corresponding encoded qubits in different blocks. To actually get universal computation, we 
will want to perform operations between qubits encoded in the same block. To do this, we 
need a few more tools, which will be presented in the next section. I will also consider a few 
more examples where we have tools beyond the ones available for any code. 

VII. DISTANCE 2 CODES 

There is a large class of distance 2 codes with a very simple form. The stabilizer for these 
codes has just two generators, one a product of all As and one a product of all Zs. The 
total number of qubits n must be even. These codes encode n — 2 qubits, and therefore serve 
as a good model for block codes encoding multiple qubits. While these distance 2 codes 
cannot actually correct a general error, they may be useful in their own right nonetheless. 
A distance 2 code can be used for error detection [17]. If we encode our computer using 



distance 2 codes, we will not be able to fix any errors that occur, but we will know if an 
error has invalidated our calculation. A better potential use of distance 2 codes is to fix 
located errors [|TJ|. Suppose the dominant error source in our hardware comes from qubits 
leaving the normal computational space. In principle, without any coding, we can detect 
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not only that this has happened, but in which qubit it has occurred. We can then use 
this information in conjunction with a distance 2 code to correct the state, as with a usual 
quantum error- correcting code. A final possible use of distance 2 codes is to concatenate 
them to produce codes that can correct multiple errors. Since the limiting factor in the 
computational threshold for concatenated codes is the time to do error correction, this 
offers potentially a great advantage. However, there is a significant complication in this 
program, since the codes given here encode more than one qubit, which complicates the 
concatenation procedure. 

Because of the simple structure of these distance 2 codes, we can immediately see a 
number of possible fault-tolerant operations. The bitwise Hadamard rotation and the bitwise 
CNOT are both permissible. If the total number of qubits is a multiple of 4, the P gate 
and the other single qubit operations are allowed, as well. What is less clear is how these 
various operations affect the encoded data. 

The Xi operators for these codes are XiX i+i , where % runs from 1 to n — 2. The Zi 
operators are Z i+ iZ n . Therefore, swapping the (i + l)th qubit with the (j + l)th qubit will 
swap the ith encoded qubit with the jth. encoded qubit. Swapping two qubits in a block is 
not a transversal operation, but if performed carefully, it can still be done fault-tolerantly 
One advantage of the swap operation is that any errors in one qubit will not propagate to the 
other, since they are swapped as well. However, applying the swap directly to the two qubits 
allows the possibility of an error in the swap gate itself producing errors in both qubits. We 
can circumvent this by introducing a third ancilla qubit. Suppose we wish to swap A and B, 
which are in spots 1 and 2, using ancilla C, in spot 3. First swap the qubits in spots 1 and 
3, then 1 and 2, and finally 2 and 3. Then A ends up in spot 2, B ends up in spot 1, and 
C ends up in spot 3, but A and B have never interacted directly. We would need two swap 
gates to go wrong in order to introduce errors to both A and B. Note that while the state 
C docs not matter, it should not be something important, since it is exposed to error from 
all three swap gates. Also note that we should perform error correction before interacting 
this block with another block, since errors could then spread between corresponding qubits, 
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which have changed. 

The action of the CNOT is simple. As for other CSS codes, it just produces a CNOT 
between all of the encoded qubits in the first block with the corresponding encoded qubits 
in the second block. The Hadamard rotation converts Xi to ZiZi + i, which is equivalent 



(via multiplication by M2) to Z 2 . . . ZiZi + 2 ■ ■ ■ Z n . This is equal to Z\ . . . Zi_\Zi+\ . . . Z n - 



2- 



Similarly, Zi becomes X\ . . . Xi-iXi+i . . . X ra _2. For instance, for the smallest case, n = 4, 

Xi — > Z2 

Z;^X~ 2 (47) 

X2 — > Z\ 

Z2 — *• X\. 

The Hadamard rotation for n = 4 performs a Hadamard rotation on each encoded qubit 
and simultaneously switches them. For larger n, it performs the Hadamard rotation on each 
qubit, and performs a variation of the class of codes discussed in section [VT| . 
For n = 4, the P gate acts as follows: 



X\ — > —Y\i2 — —X1Z2 



X 2 -> -Y 1 Y 3 = -X 2 Z, (48) 

Z~i^Z~x 

Z2 — > Z2- 

A consideration of two-qubit gates allows us to identify this as a variant of the conditional 
sign gate. Specifically, this gate gives a sign of —1 unless both qubits are |0). 

When we allow measurement, a trick becomes available that is useful for any multiple- 
qubit block code. Given one data qubit, prepare a second ancilla qubit in the state |0) + |1), 
then apply a CNOT from the second qubit to the first qubit and measure Z for the first 
qubit. The initial stabilizer is I <S> X; after the CNOT it is X <g> X. Therefore the full 
operation takes X <g> / to / <g> X and Z <g> / to Z <g> Z. We can discard the first qubit and the 
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second qubit is in the initial data state. However, if we prepare the ancilla in the state |0), 
then apply a CNOT, the original state is unaffected. Therefore, by preparing a block with 
all but the jth encoded qubit in the state |0), and with the jth encoded qubit in the state 
|0) + |1), then applying a CNOT from the new block to a data block and measuring the jth 
encoded qubit in the data block, we can switch the jth encoded qubit out of the data block 
and into the new, otherwise empty block. 

This trick enables us to perform arbitrary operations on qubits from the same block for 
the distance 2 codes. We switch the qubits of interest into blocks of their own, use swap 
operations to move them into corresponding spots, then perform whole block operations 
to interact them. Then we can swap them back and switch them back into place in their 
original blocks. 

The step that is missing for arbitrary stabilizer codes is the ability to move individual 
encoded qubits to different places within a block. Since the gate in section [VI] gives us 
a block CNOT, we can perform the switching operation into an empty block. By using 
switching and whole block operations, we can perform an arbitrary one-qubit operation on 
any single encoded qubit within a block. The only remaining operation necessary is the 
ability to swap an encoded qubit from the ith place to the jth place. We can do this using 
quantum teleportation. All that is required is an otherwise empty block with the ith and jth 
encoded qubits in the entangled state |00) + |11). Then we need only perform single-qubit 
operations and a CNOT between the qubits in the ith places, both of which we can do. 
To prepare the entangled state, we simply start with the +l-eigenstate of Zi and Zj, then 



measure the eigenvalue of X t Xj (and correct if the result is —1). This is just an operator in 



Q, so we know how to do this. The state stays in an eigenvector of ZiZj, which commutes 



with XiXj, so the result will be the desired encoded Bell state. We can then teleport the zth 
qubit in one otherwise empty block to the jth qubit in the block originally containing the 
entangled state. This was all we needed to allow universal computation on any stabilizer 
code. 
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TABLE III. The stabilizer and encoded Xs and Zs for the eight-qubit code. 
VIII. THE 8 QUBIT CODE 



There is a code correcting one error encoding 3 qubits in 8 qubits P,|5|,[T9|j . The stabilizer 
is given in table [TIT]. There are no transversal operations that leave this stabilizer fixed. 
However, when we allow swaps between the constituent qubits, a number of possibilities 
become available. 

One possible operation is to swap the first four qubits with the second four qubits. This 
leaves Mi, M 2 , and M4 unchanged. M3 becomes instead M1M2M3, and M5 becomes M1M5. 
On the encoded qubits, this induces the transformation 

X <8> I® I ->X ®I ® Z 
I ®X ®I -»•/ ®X®I 

I®I®X -> Z®I®X (49) 

Z ® I ®I -»• Z ®I ®I 

I®Z®I^I®Z®I 
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I ®I ®Z -> I ®I ®Z. 

This is just a conditional sign on the first and third qubits, with the second encoded qubit 
unaffected. Through single-qubit transformations, we can convert this to a controlled NOT, 
and using this perform a swap between the first and third encoded positions. 

Another operation is to swap qubits one and two with three and four and qubits five 
and six with seven and eight. This leaves Mi, M2, and M3 unchanged, and converts M4 to 
M2M4 and M5 to M1M5. On the encoded qubits, it induces the transformation 

X®I®I->X®Z®Z 

I ®X ® I ->■ z ®x ® z 

I®I®X -»■ Z®Z®X (50) 

Z ® I ® / ->■ Z <g) I <g> I 
I®Z®I->I®Z®I 
I ® I ® Z -> I ® I ® z. 

We could also switch the odd numbered qubits with the even numbered qubits. That leaves 
Mi and M 2 unchanged, while turning M 3 into M t M 3 , M 4 into M t M 4 , and M 5 into M x M 2 M b . 
On the encoded qubits it induces 

X ® J ® J -> X® I®Z 

I ®X ® I ^ I ®x ® z 

I®I®X ^ Z®Z®X (51) 

Z ® I ® I -> Z ®I ® I 

I®Z®I^I®Z®I 
I®I®Z^I®I®Z. 

This is just a conditional sign between the first and third places followed by a conditional 
sign between the second and third places. Combined with the first operation, it gives us a 
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conditional sign between the second and third places, which we can again convert to a swap 
between the second and third encoded positions. This allows us to swap any two encoded 
qubits in the block, which is sufficient to give us universal computation. 

In this case, the symmetries of the code naturally became allowed transformations of the 
stabilizer. This is likely to hold true in many other cases as well. As with the five-qubit 
code, we also had available a universal protocol for swapping encoded qubits, but multiple 
methods again allow us more freedom in choosing efficient methods. 

IX. SUMMARY AND DISCUSSION 

I have presented a general theory for understanding when it is possible to apply a given 
operation transversally to a given quantum error-correcting code, and for understanding the 
results of making a measurement on a stabilizer code. These results clarify the advantages 
of the doubly-even self-dual CSS codes used by Shor 0. They also provide protocols for 
performing universal computation on any stabilizer code. In many cases, the protocols 
described here call for a number of steps to perform most simple operations, so more efficient 
protocols for specific codes are desirable, and I expect the methods described in this paper 
will be quite helpful when searching for these protocols. 

Efficient use of space is also important. Existing methods of fault-tolerant computation 
use space very inefficiently, and being able to use more efficient codes (such as those encoding 
multiple qubits in a block) could be very helpful in reducing the space requirements. 
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by the Army Research Office. I would like to thank John Preskill, Manny Knill, Richard 
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